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Abstract 

The minimal supersymmetric standard model with soft breaking has a large landscape of su- 
persymmetric particle mass hierarchies. This number is reduced significantly in well-motivated 
scenarios such as minimal supergravity and alternatives. We carry out an analysis of the landscape 
for the first four lightest particles and identify at least 16 mass patterns, and provide benchmarks 
for each. We study the signature space for the patterns at the CERN Large Hadron Collider by 
analyzing the lepton -|- (jet > 2) -|- missing Pt signals with 0, 1, 2 and 3 leptons. Correlations in 
missing Pt are also analyzed. It is found that even with 10 fb~^ of data a significant discrimination 
among patterns emerges. 

PACS numbers: 12.60.Jv,04.65.+e,11.30.Pb,14.80.Ly 
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Introduction: The search for supersymmetric particles (sparticles) is one of the major 
goals of the current experiments at the Fermilab Tevatron collider and at the CERN Large 
Hadron Collider (LHC) which will come on line in the very near future. Central to the 
discovery of sparticles is the way their masses align in a hierarchical mass pattern as such 
alignment has strong influences on the ability of experiments to probe signatures emerging 
from sparticle production at colliders. There are 32 supersymmetric masses in the minimal 
supersymmetric standard model (MSSM) including the Higgs bosons. With the sum rule 
constraints in the gaugino and Higgs sectors but without any phenomenological constraints 
the number of mass parameters is 27. Using Stirling's formula [n! ~ \/27m(n/e)"] one finds 
at least 0(10^'^) possibilities for hierarchical patterns. Although this number does not rise 
to the level of O(10^°°°) as in string landscapes, it is still an impressive number and can 
be construed as a mini landscape in this context. The number of possibilities is drastically 
reduced in the minimal supergravity grand unified model mSUGRA although no 

classification has ever been made and the precise number of possibilities is not known. In 
this Letter we undertake this cartography. To keep the analysis under control we focus on 
mass hierarchies for the first four sparticles excluding the lightest Higgs. Such a cartography 
is important for devising strategies for analyzing data from the LHC. 

Our analysis, based on a large scan of the mSUGRA parameter space, reveals at least 16 
hierarchical patterns for the first four sparticles, of which a significant number do not appear 
in conventional benchmarks such as Snowmass points (these are a representative sample 
of points in the MSSM parameter space called benchmarks which are used for theoretical 
jredictions of supersymmetry at present and future colliders) and Post WMAP benchmarks 
^ (which is an updated version of previous benchmarks taking account of the relic density 
constraints from data from the Wilkinson Microwave Anisotropy Probe (WMAP)). The 

n n u 

analysis is carried out in the mSUGRA model (for recent works see |5|, |6|, 17J). mSUGRA 
is an effective theory below the Grand Unification scale Mq and its parameter space is 
defined by the parameters (mo, mi/2, Ao, tan/3, sign/x) where (mo,mi/2) are the universal 
scalar and gaugino masses, Aq is the universal trilinear coupling, and tan/5 is the ratio of 
the vacuum expectation values (VEVs) of the two neutral Higgs fields in MSSM, and sign/^ 
is the sign of /i, where fi is the Higgs mixing parameter. In the analysis we impose the 
relic density constraints on the abundance of the lightest neutralino (x?) consistent with 
the WMAP: 0.0855 < Q^oh'^ < 0.1189 jsl, and other experimental constraints as follows: 
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FIG. 1: The surviving hierarchical patterns in the mSUGRA landscape with constraints as dis- 
cussed in the text. In mSP14,15,16 the Xi the Higgs bosons {A, H) are sometimes seen to 
switch. 

2.83 X 10"^ < BR(b 37) < 4.63 x lO""^ (where the SUSY loop contributions to this 
decay can be comparable to the corrections in the Standard Model (SM)), nih > 100 GeV, 



m 



~± > 104.5 GeV, m-^^ > 101.5 GeV, m:?^ > 98.8 GeV, where h,xf,ti,fi are the lightest 
Higgs boson, chargino, stop and stau respectively. _^ 

For the calculations of the relic density of Xi we use MicrOMEGAs version 2.0.1 | |9|] with 
sparticle and Higgs masses calculated using the package SUSPECT 2.3 of Ref. [lO|]. We 



have investigated other softwares jll. [l2. 
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15|, 



30 



and find significant agreement 



among them. For each mSUGRA model point that survives the constraints mentioned above, 
we compute the sparticle spectrum, mixing angles, etc. The resultant SUSY Les Houches 
Accord (SLHA) [18] file is called by the PGS4 Olympics main Fortran file, which simulates 



19| for detailed discussion). For the 



the LHC detector effects and event reconstruction (see 
computation of SUSY production cross sections and branching fractions we employ PYTHIA 
6.4.11 20(], which simulates hadronization and event generation from the fundamental SUSY 



Lagrangian. SUSY cross section and events were generated using all 32 sparticles in this 
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1785 mSUGRA models out of 2 million point scan with Monte Carlo simulation 
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FIG. 2: The plots exhibit the 16 mass patterns (see Fig.(l)) in tan/3 — Ao/niQ and ttiq — 
planes among the 1785 surviving points out of a random scan using 2 x 10® models with flat priors 
in the ranges indicated in the left figure panel consistent with the constraints discussed in the text. 



analysis. We have cross checked sample points with PROSPINO 2l|] which computes the 
next-to-leading order cross sections for the production of supersymmetric particles at hadron 
colliders, and TAUOLA [22] is called by PGS4 for the calculation of tau branching fractions. 
With PGS4 we use the Level 1 (LI) triggers based on the Compact Muon Solenoid detector 
(CMS) specifications [23j and the LHC detector card. 

Landscape of mass hierarchies: We focus on the patterns for the four lightest particles 
(discounting the lightest Higgs) as they would to a great degree influence the discovery of 
SUSY while keeping the size of the landscape in check. We have carried out a mapping of 
the mSUGRA parameter space with 2 x 10® models with /i > and the parameter range 
given in Fig. ([2]) via Monte Carlo scan with flat priors under the constraint of radiative 
breaking of the electroweak symmetry. 55 hierarchical patterns for the first four particles 
were seen which are reduced to 16 with relic density and collider constraints. We label these 
as mSUGRA pattern 1 (mSPl) through mSUGRA pattern 16 (mSP16) as shown in Fig.([2]). 
The frequency of their occurrence is exhibited in Fig.([l]). A significant set of the mSPl 
points lie in the region |y4o/mo| < 1 and correspond to the Hyperbolic Branch/ Focus Point 
(HB/FP) region 2J]. In Tabled we also give illustrative, mass pattern motivated, benchmark 
points, for each of the mass patterns mSPl-mSP16 chosen with a moderate to light SUSY 
scale {Q = ^Jfrv^Jni^) and to show the effects of scanning over the parameter space. Most 
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TABLE I: Benchmarks using SUSPECT 2.3 with one point for each mass pattern mSPl-mSP16. 
Also given are the neutrahno LSP (Lightest SUSY (R parity odd) Particle), and the Lightest 
Charged Particle (LCP) masses. We take /i > 0, mb^{mb) = 4.23 GeV, a^^(Mz) = .1172, 
and mt(pole) = 170.9 GeV. At least five LCP from these benchmarks will be accessible at the 
International Linear Collider (ILC). 
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of the mSP patterns do not appear in previous works. Thus all the Snowmass mSUGRA 
points (labeled SPS) jsl are only of types mSP(l,3,5,7) as follows: (SPSla, SPSlb, SPS5) 

^ mSP7, SPS2 ^ mSPl, SPSS ^ mSP5, and (SPS4, SPS6) ^ mSP3. Regarding the 

j 

Post-WMAP points |4| one has the following mapping (A', B', C, D', G', H', J', M') ^ mSP5, 
(I', L') mSP7, E' mSPl, K' — > mSP6. We observe more patterns owing to the sampling 
of larger regions of the parameter space especially in A^. 
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FIG. 3: Leptons + n jets {n > 2) signatures originating from 16 patterns with Pt and rapidity rj 
cuts as shown where (a) single-lepton vs no-lepton; (b) single-lepton plus single-r vs dilepton. the 
shaded regions in (a)-(b) are due to SM backgrounds it, 66, Dijets, Drell-Yan, and Z, W production. 



The discovery hmits are max.{b\J Nsm ^ 10}, where Nsm is the Standard Model background. 
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FIG. 4: The number of r^r~ events vs (a) trilepton events, and vs (b) P™*^ cut. In (b) the post 
trigger level cuts are as in Figs. ([3]) and (4a) except for the PJp"^^" cut. Each curve is for one point 
for each mSP, chosen such that their PGS4 data files are nearly of the same (5MB) size. 



Mass hierarchies and their signatures: We have carried out analyses of the lepton + 
jet signals to determine if such signals can discriminate among the patterns. We have 
analyzed 0, 1, 2 and 3 leptonic events (lepton = L, r; L = e, yu) with and without jets. 
Our Standard Model backgrounds are checked against two CMS notes [25| and our results 



lie within the error bars of these analyses. Some analyses of this type exist 



27 



28|, 



29| 
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but not in the context of mass hierarchies analyzed here. We give now the details of the 
analysis. In Fig.Q, we use 902 sample models from Fig. ([2]) (corresponding to a scan of 10^ 
points as it contains all the essential features of the larger scan and is visually clearer) to 
generate SUSY signals through PGS4 using 10 fb^^ of integrated luminosity. In Fig.([3k) an 
analysis of the IL + (jets > 2) vs OL + (jets > 2) events is given along with the discovery 
limit for each of these model points. Here one is beginning to see discrimination among 
patterns. Specifically the signals of mSPll-mSP13 lie significantly lower than the others. 
In Fig.([3]D) we give an analysis of Ir + IL + (jets > 2) vs 2L + (jets > 2). We see that 
relative to Fig.(l3tL) there is now more dispersion among the patterns and furthermore the 
SM background is also significantly smaller in this case. Finally in Fig.(j3^) we have a plot 
of OS(Opposite Sign)2r + (jets > 2) vs 3L + (jets > 2). Here we see that the dispersion 
among models is even larger and the SM background is even smaller. Thus we can see 
visually the effect of the SM background shrinking, the signal relative to the background 
getting stronger and the discrimination among the patterns increasing as we move from 
Fig.(l3ti) to Fig.(jlti). As expected the trileptonic signal [30| is the strongest with the least 
background. The missing Pt cuts can also help discriminate among models. In Fig. (jib) we 
give an analysis of OS 2r + (jets > 2) vs missing Pt cut. One may note the very significant 
dispersion among the patterns with missing Pt cut. 

Conclusion: The discovery of supersymmetry is one of the major goals of the current 
effort at the Tevatron and in experiments with the CMS and the ATLAS detectors in the 
very near future at the LHC. When sparticles are produced the signatures of their production 
will be determined by their hierarchical mass patterns. In this Letter we have investigated 
hierarchical mass patterns for the four lightest sparticles within one of the leading candidate 
theories - the mSUGRA model. The analysis shows 16 such patterns consistent with all the 
current experimental constraints and most of these patterns are new and not discussed in 
the previous literature. We also carried out an analysis of signatures of these patterns which 
include (0,1,2,3) leptons + n jets (n > 2) signals and also missing transverse momentum 
{Pt) signals. We conclude that even with 10 fb^^ of luminosity significant dispersion can 
be seen among many mSP signatures and thus discrimination among patterns is possible. 
Similar analyses are desirable for other plausible soft breaking scenarios beyond mSUGRA. 

Acknowledgements: We thank Brent Nelson and George Alverson for discussions. This 
work is supported in part by NSF grant PHY-0456568. 



7 



[1] A. H. Chamseddine, R. Arnowitt and P. Nath, Phys. Rev. Lett. 49, 970 (1982); P. Nath, 
R. Arnowitt and A. H. Chamseddine, Nucl. Phys. B 227, 121 (1983). 

[2] R. Barbieri, S. Ferrara and C. A. Savoy, Phys. Lett. B 119, 343 (1982); L. J. Hall, J. D. Lykken 
and S. Weinberg, Phys. Rev. D 27 (1983) 2359. 



[3] B. C. Allanach et al, |arXiv:hep-ph/0202233| . 

[4] M. Battagha, A. De Roeck, J. R. Ellis, F. Gianotti, K. A. Olive and L. Rape, Eur. Phys. J. C 
33, 273 (2004). 

[5] A. Djouadi, M. Drees and J. L. Kneur, JHEP 0603, 033 (2006); U. Chattopadhyay, D. Das, 
A. Datta and S. Poddar, Phys. Rev. D 76, 055008 (2007); M. E. Gomez, T. Ibrahim, P. Nath 
and S. Skadhauge, Phys. Rev. D 72, 095008 (2005). 

[6] R. Arnowitt et al, Phys. Lett. B 649, 73 (2007); V. Khotilovich, R. Arnowitt, B. Dutta and 
T. Kamon, ibid, 618, 182 (2005). 

[7] J. R. Elhs, S. Heinemeyer, K. A. Olive and G. Weiglein, JHEP 0605, 005 (2006). 



[8] D. N. Spergel et a^., !! astro-ph/0603449[ 

[9] G. Belanger, F. Boudjema, A. Pukhov and A. Semenov, Comput. Phys. Commun. 176, 367 
(2007). 

[10] A. Djouadi, J. L. Kneur and G. Moultaka, Comput. Phys. Commun. 176, 426 (2007). 



[11] F. E. Paige, S. D. Protopopescu, H. Baer and X. Tata, |arXiv:hep-ph/0312045 

[12] W. Porod, Comput. Phys. Commun. 153, 275 (2003). 

[13] B. C. Allanach, Comput. Phys. Commun. 143, 305 (2002). 

[14] H. Baer, J. Ferrandis, S. Kraml and W. Porod, Phys. Rev. D 73, 015010 (2006). 

[15] G. Belanger et.al., Phys. Rev. D 72, 015003 (2005). 

[16] B. C. Allanach, A. Djouadi, J. L. Kneur, W. Porod and P. Slavich, JHEP 0409, 044 (2004). 



[17] B. C. Allanach, S. Kraml and W. Porod, JHEP 0303, 016 (2003). and |arXiv:hep-ph/0207314 
[18] P. Skands et al, JHEP 0407, 036 (2004). 

[19] |http : //www . physics . ucdavis . edu/~coiiway/resecirch/sof tware/pgs/pgs4-general . htm 

[20] T. Sjostrand, S. Mrenna, P. Skands, JHEP 0605, 026 (2006). 
[21] W. Beenakker et.al., Phys. Rev. Lett. 83, 3780 (1999). 

[22] S. Jadach, Z. Was, R. Decker and J. H. Kuhn, Comput. Phys. Commun. 76, 361 (1993). 



8 



[23] CMS Collaboration, CERN/LHCC 2006-001 (2006). 

[24] K. L. Chan, U. Chattopadhyay and P. Nath, Phys. Rev. D 58, 096004 (1998); J. L. Feng, 
K. T. Matchev and T. Moroi, Phys. Rev. Lett. 84, 2322 (2000); H. Baer, C. Balazs, A. Belyaev, 
T. Krupovnickas and X. Tata, JHEP 0306, 054 (2003). 

[25] D. J. Mangeol, U. Goerlach, CMS NOTE 2006/096, M. ChiorboH, M. Galanti, A. Tricomi, 
CMS NOTE 2006/133 

[26] N. Arkani-Hamed, G. L. Kane, J. Thaler and L. T. Wang, JHEP 0608, 070 (2006). 

[27] J. P. Conlon et.al., JHEP 0708, 061 (2007). 

[28] G. L. Kane, P. Kumar and J. Shao, J. Phys. G 34, 1993 (2007) 

[29] H. Baer, V. Barger, G. Shaughnessy, H. Summy and L. t. Wang, Phys. Rev. D 75, 095010 
(2007). 

[30] P. Nath and R. Arnowitt, Mod. Phys. Lett. A 2, 331 (1987); R. Arnowitt, R. M. Barnett, 
P. Nath and F. Paige, Int. J. Mod. Phys. A 2, 1113 (1987); H. Baer, C. h. Chen, F. Paige and 
X. Tata, Phys. Rev. D 50, 4508 (1994); V. D. Barger, C. Kao and T. j. Li, Phys. Lett. B 433, 
328 (1998). 



9 



